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g> ' Abstract 

We study the one-loop contributions of matter and radiation to the grav- 

in ' 

| itational polarization tensor at finite temperatures. Using the analytically 

continued imaginary-time formalism, the contribution of matter is explicitly 
given to next-to-leading (T 2 ) order. We obtain an exact form for the con- 
| tribution of radiation fields, expressed in terms of generalized Riemann zeta 

On" 

<^ ' functions. A general expression is derived for the physical polarization tensor, 

which is independent of the parametrization of graviton fields. We investi- 
gate the effective thermal masses associated with the normal modes of the 
corresponding graviton self-energy. 



Oh 



■ I. INTRODUCTION 



Many properties of plasmas in thermal field theories can be understood from the study 
of the polarization tensor evaluated at finite temperature Jl|-§[. This tensor, which is the 
two-point correlation function, describes phenomena such as the propagation of waves and 
dumping of fields in the plasma. In thermal quantum gravity, the behavior of the polarization 
tensor is also of interest, especially in connection with cosmo logical applications. If the 
temperature T is well below the Planck scale, perturbation theory can be used to calculate 
the thermal Green functions. Thus, one obtains loop-diagrams in which the internal lines 
represent matter and radiation in thermal equilibrium, and the external lines represent the 
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gravitational fields. There has been a lot of work on hot quantum field theory in the presence 
of a gravitational field Thus far these investigations have been mainly restricted to the 

study of the hard thermal loops contributions, which are obtained in the high temperature 
limit. 

The purpose of this work is to study the behavior of the graviton polarization tensor 
at all temperatures, which might be useful in some applications. Since these calculations 
are considerably more complicated than those performed at high temperatures, we have 
restricted for definiteness to work to one-loop order with thermal bosonic fields, which may 
be of spin or 1. The method we use is that of reference 0, where the Green functions 
are related to a momentum integral of the forward scattering amplitude of thermal particles 
in a gravitational field. Then, the temperature-dependent part of the graviton polarization 
tensor can be written at all temperatures in the form: 

U ^,aP / fc ) = _J_ f A 1_ Fmaf > / k s ,j A , 

K } (2n) 3 J 2Qexp(Q/T)-l w ' K J 

Here g M = (Q, q) represents the on-shell momenta of a thermal particle with mass m and 
energy Q = \f\q\ 2 + m 2 . F^ v ' al3 (q, k) is the forward scattering amplitude, summed over 
the polarizations of thermal particles, which is a covariant function of q and the external 
momenta k. This temperature-independent amplitude is weighted in (|1 . 1|) by the Bose dis- 
tribution factor. Because of the angular integrations, H^"' 01 * is no longer a Lorentz covariant 
function. It depends on the time-like vector u^, representing the local rest frame of the 
plasma. For simplicity, we work in the comoving coordinate system where = 5q. The 
above method simplifies very much the calculations in the present case. 

In Sec.[TI] we consider the contribution of matter particles described by the scalar field <f), 
coupled to a gravitational field. The coupling characterized by the term ^R<p 2 is included, 
where £ is a numerical factor and R denotes the Ricci scalar. We verify that IP 1 ^ satis- 
fies the Ward identity which reflects the invariance of the action under general coordinate 
transformations. We obtain a general expression for the leading (T 4 ) and next-to-leading 
(T 2 ) contributions to the graviton polarization tensor. The special case when £ = —1/6 and 
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m = is of particular interest, since then the scalar action is also invariant under conformal 
transformations JTIJ. Due to this invariance, IP 1 ^ satisfies in this case a Weyl identity 
which is explicitly verified. 

In Sec.[TTT| we discuss the coupling of radiation fields which may be photons or gluons, to a 
gravitational field. This coupling is also invariant under general coordinate transformations 
as well as under conformal transformations. We remark that the thermal contributions 
associated with internal gauge fields represent gauge-invariant quantities. The Ward and 
Weyl identities determine uniquely the (T 4 ) contributions, which are the same for all thermal 
particles, apart from numerical factors which count the number of degrees of freedom. Using 
general properties of the forward scattering amplitude, we show that all other contributions 
can be expressed in terms of just 2 parameters which are not fixed by the Ward and Weyl 
identities. Rather, these parameters depend specifically on the nature of thermal particles. 



In Sec.|Iy| we obtain a closed form expression for the contributions of thermal radiation 
fields to the graviton polarization tensor. We show that these can be expressed in terms 
of generalized Riemann zeta functions ((—n,t) ]TT[ for natural values of n, t being a ratio 
of external momenta and the temperature. In the high temperature limit, this expression 
yields a series of decreasing powers in the temperature, which includes leading (T 4 ) and- 
next-to-leading (T 2 ) contributions. Some technical aspects which arise in the calculations 
are discussed in the Appendices. 

In Sec.|V] we analyze the dependence on the parametrization of the graviton fields, of the 
one-particle irreducible (1PI) contributions to the graviton polarization tensor. This behav- 
ior occurs generally because of the non-vanishing of the thermal graviton 1-point function. 
We show that the physical polarization tensor, identified with the graviton self-energy, is de- 
scribed by a traceless function which includes contributions from thermal 1-point functions. 
A general expression for the physical self-energy at finite temperature is derived, which is 
independent of the graviton parametrization. 

In Sec.|VT| we discuss the effective graviton propagator, obtained by iterative insertions 
in the free propagator of the physical self-energy. We analyze, in the static limit, the 



corresponding poles which describe three normal modes of dynamical screening. While one 
of the modes remains unshielded, a no n- vanishing screening mass m 2 s = 32nGp/3 appears in 
the spatially transverse one, where p is the thermal energy density. The spatially longitudinal 
mode is characterized by an imaginary mass m 2 j = —32nGp, similar to the classical Jeans 
mass, indicating an instability of thermal quantum gravity. 



II. MATTER CONTRIBUTIONS TO THE POLARIZATION TENSOR 

We consider here thermal matter represented by scalar particles of mass m, coupled to 
the gravitational field via the Lagrangian: 

C(x) = -\sj-g{x) [gTd^drf - (m 2 + £R) 4> 2 } , (2.1) 

and expand the metric tensor g^ u in terms of the deviation from the Minkowski metric r]^: 

g^u = Vvv + Khftv, (2.2) 

where k = \f32irG. In order to derive the one-particle irreducible (1PI) contributions to 
the thermal graviton 2-point function, we consider the Feynman graphs shown in Fig.|l|. 
According to Eq. (|1 . 1|) , these can be expressed in terms of the forward scattering amplitude 
of on-shell scalar particles, as indicated in Fig.[| The corresponding contributions to the 
amplitude can be expressed in terms of a basis of 14 independents tensors Tj iu ' al3 (q, k), which 
are symmetric under the interchanges (p <-> v\ [a <->• (3) and (p, v) <-» (a,/3). These tensors 
are covariant functions of q and k, being polynomials of maximum degree 4 in the momenta. 
They can be obtained from Table H, replacing the vectors (X, Y) by the pair (q, k). With help 
of the Feynman rules given in Appendix 0, it is straightforward to obtain for the forward 
scattering amplitude the expression: 
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If" °* + (i £ A; 4 + e 2 fc 4 - ^ 2 fc " 9 + i(* • qf) 

(Ik 2 + £ k 2 ) Tp a( " + hrp aP + (h k 2 + ^k-q\ Tp aP - 

S 9 ^ 10 2 ^ 

(~£ * 2 + e 2 ^ 2 ) T/T 0/3 + (~£ A: 2 - i* • g) Tff* a/J 
+ (A; <-> -A;) . 

(2.3) 

To obtain the leading (T 4 ) and the next to leading (T 2 ) contributions to the polarization 



tensor, we need to expand the energy Q = y\q\ 2 + m 2 in powers of (m 2 / |g| 2 ), as well as the 
Feynman denominators: 

11 k 2 k 4 k 6 

k 2 + 2k ■ q ~ 2k ■ q ~ (2k ■ qf + (2k ■ qf ~ (2k ■ qf + " ' ' ^' ' 

The T 4 contributions come from terms in the forward amplitude ( |2.3|) which are homoge- 
neous functions of q of degree 2. These are given by: 

-T^ (q, k) - -^TjTr^ (q, k) + -Llf"** (k, q) = 
(k-q) k-q 

k 2 

- (rj^q^q 01 + r/^q^q 13 + rj^q v q a + rj^q^qn ^q f "q u q a q 13 ■ (2.5) 

^ ' (k-q) 

(q^q v q a k p + q»q v k a q p + q t *k v q a q p + k»q v q a q p ) 
k ■ q v 1 

Note that terms involving the parameter £ do not contribute to Q2.5|) . These contribute only 

to next-to- leading (T 2 ) order, which result from terms of degree zero in q in the forward 

amplitude. In order to find these contributions, we perform the \q\ integration in ( |1 . 1|) using 

the formulas: 

\%m _^ 2 



o [exp(|g1/T)-l] 6 
qfd\q\ vr 4 T 4 vr 2 m 2 T 2 



Q [exp (Q/T) - 1] 15 4 



+ •••. (2.7) 
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The angular integrals can be done using the methods described in 0. The result can be 
expressed in terms of the basis of 14 tensors Tj iu,al3 (u, K), obtained from Table |, where we 
replace the pair (X,Y) by (u,K). Here = 5q and 

k>* ( k 



(r, k) . 



(2.8) 



Then the 1PI contributions to the polarization tensor can be written up to the next to 
leading order in the form: 
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H^> <# (k, m, = £ II* (r, K, £) If' aP («, K) 



(2.9) 



i=i 



where: 
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tt 2 T 4 h (r, K) + k T n t (r, K, f) + m 2 T 2 Sl (r, K, f) 



+ 



(2.10) 



The explicit form of the dimensionless functions Zj(r, K), rii(r, K,£) and Sj(r, K, £) are given 
in Appendix ||. These exhibit, apart from a logarithmic dependence in r, a polynomial 
behavior in K of maximum degree 10. The coefficients U (r, K) which contribute to the 
leading (T 4 ) order, have been obtained previously || and are included here for completeness. 

As a consequence of the invariance of the theory under general coordinate transforma- 
tions, the 1PI graviton 2-point function satisfies the Ward identity: 



-k v IF"- aP (k) = F T af3 - k a (r aa if* + T Pa r) a ") . 
Here r a/3 denotes the thermal graviton 1-point function, which is given by: 



(2.11) 



7T 2 T 4 K 



m 2 T 2 K 



(2.12) 



With the help of the expression given by Eqs. ( |2.9|) and ( p.lO| ), the Ward identity (|2.11|) can 
be explicitly verified to this order. It is well known [|K]] that in the conformally coupled case, 
when £ = —1/6 and m = 0, the action is also invariant under conformal transformations 
given by: 



(x) = fl 2 (x) g^, 



(2.13) 



(x) = JT 1 (x) 4> (x) . (2.14) 

In consequence of this invariance, the 1PI graviton 2-point function will also satisfy the Weyl 
identity @,|]: 

-W v >%{k) = -T» v . (2.15) 

K 

This identity is explicitly verified by our expression for IP^"^ (k) [Eq. ( |2.9|) 1 and for T^ v 
[Eq. fl2TP2D 1 evaluated at £ = -1/6 with m = 0. 

III. RADIATION FIELDS CONTRIBUTIONS TO THE POLARIZATION 

TENSOR 

In this section we analyze the contributions of spin 1 gauge fields, which may be photons 
or gluons. Since for our purpose the self-interactions of the Yang-Mills particles can be 
neglected, there is no loss of generality in considering only the contribution of an Abelian 
field A^" . For non- Abelian fields the contributions are the same, up to an overall color factor. 
The coupling of the gauge field A 11 is described by the Lagrangian: 

C A = --J-g(x)g^g^ (d^A a - d a AJ {d u A p - d p A v ) . (3.1) 

It is convenient for computational purposes to fix the gauge by choosing: 

= (V,^) (V^) , (3.2) 

where is the covariant derivative. The corresponding Faddeev Popov Lagrangian is given 
by: 

Up = gT\l-9 (x) (d.x) (duX) , (3.3) 

where \ and \ are the ghost fields. The form of the above interactions is such that the 
theory is invariant under local coordinate transformations, as well as under conformal trans- 
formations given by: 
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(x) = (x) ; (x) = Q 2 (x) g^ v . (3.4) 

As we have seen, these invariances ensure the 1PI graviton 2-point function to satisfy the 
Ward and Weyl identities given respectively by Eqs. Q2-ll|) and (|2.15| ). Here T a ^ is obtained 
multiplying (|2.12|) by a factor 2 and setting m = 0. 

With help of the Feynman rules listed in Appendix |A], we can evaluate the 1PI graphs 
contributing to rP^ a/3 , which are shown in Fig.|3[ The diagrams contributing to the cor- 
responding forward amplitude are represented in Fig.[|. It is important to note that the 
thermal contributions from internal gauge fields represent gauge-independent quantities. 
We have verified this independence explicitly, performing all computations in the general 
class of covariant gauges defined by ( |3T2"|) . The dependence on the gauge parameter a cancels 
in the final expression of the forward scattering amplitude, which is given by: 

-\{k-qf (q, k) - k \t a " (q, k) + (?, k) 

+l T r a(3 (q, k) + l -T^ a " (q, k) - ^Tfr? (q, k)\ 
+ (k <-> —A;) 

(3.5) 



At this point, it is interesting to compare ( |3.5|) with the amplitude corresponding to the 
scalar case [Eq. 02.31) 1, evaluated for m = and £ = —1/6. We see that in both amplitudes, 
the coefficients of the tensors Tj*"' (i = 2, 3, 7) which contribute in the high temperature 
limit [cf. Eq. (|2.5| )] are the same, up to a factor of 2 which counts the degrees of freedom 
of a physical gauge particle. On the other hand, all other coefficients seen to be different in 
general. 

In order to understand this behavior we consider now the consequences of the Ward 
( 2.111 ) and Weyl Q2.15 ) identities on the structure of the forward scattering amplitudes. To 



this end, we use the following representation of the graviton one-point function: 

1 " (27T) 3 J Q exp (Q/T) - 1 q q ■ {6 - b) 
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Then, we find from Eq. fll.lj ) that the Ward and Weyl identities ensure the forward scattering 
amplitudes to obey respectively the relations: 

J-k u F^' aP (q, k) = Wq a q f) - k ■ q ( q a r]^ + 9 V") , ( 3 - 7 ) 

L F ^ a a (q,k) = -2q^ (3.8) 

We will now investigate the constraints imposed by the relations ( |3.7|) and ( |3.8| ) on the 
general form of the amplitudes F^ v,OL ^ . Since these are Lorentz covariant functions of q and 
k, they can be expressed in terms of the tensor basis Tf v ' a/3 (q, k) as follows: 
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F^oP (q, k) = n 2 Yl F i k ■ ?) (?. k ) » ( 3 - 9 ) 

i=X 

where Fi are invariant functions of k 2 and k ■ q. Inserting ( |3.9| ) into the Ward identity ( |3.7| ) 
and identifying the coefficients of the independents tensor structures yields 10 relations 
among the F^. Similarly the use of the Weyl identity ( |3.8| ) gives 4 more relations. However, 
not all of these relations are independent, so that we can express 11 functions Fj in terms 
of the remaining 3 as follows: 

F x 
F 2 
F 3 

F, 
F 5 

F 6 
F 7 

F 8 
F 9 



3k 4 

—F 12 + k ■ qF 14 


(3.10a) 


1 k 4 k 2 

2 2k ■ q k-q 


(3.10b) 


k 2 2k 6 k 4 
2(k-q) 2 (k-q) 3 11 ' (k-q) 3 14 


(3.10c) 


- k -F X2 - (k ■ q)F u 


(3.10d) 


k2 F 
(k-q) Fu 


(3.10e) 


k 2 

— Fu — Fxi 


(3.10f) 


1 2k 4 k 2 
2k -q (k-q) 2 " (k-q) 2 14 


(3.10g) 


k-q 3 k 2 
2 Fn 4 Fl2 


(3.10h) 


2k 2 2 

Fxx 1 , Fx4 
k-q k-q 


(3.10i) 
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k 2 

F13 = —F 12 (3.10k) 

We see that the invariance of the theory under local coordinate and conformal transfor- 
mations does not fix the functions Fn, F12 and F\±. Further constraints are provided by 
the property of the forward scattering amplitude of being a function with dimension of 
(momenta) 2 , which is even under (k < > —k). Furthermore, to one loop order in perturba- 
tion theory this amplitude can have at most one denominator involving (k 2 ± 2k ■ q). For 
instance, these general properties require the functions F 3 and F 7 to have the structure: 

F3 = C3 {Wlkk—q + W^2k— q ) ' (3 - H) 

Fl = Cl {k 2 +\k. q - k 2 -2k-q)' (3 ' 12) 

where c 3 and C7 are constants. Furthermore, it follows that F 2 must be an even function of 
k, having the structure: 

F2 = ^ ' q {wT2k— g ~ W^2k— q ) + ^ {wT2k— q + 1^2k— q ) ' (3 - 13) 

where c 2 and c' 2 are constants. Similar structures can be found for all other functions 
appearing in Eqs. ( |3.10|) . These structures yield a set of relations which must be satisfied 



identically in Eqs. ( |3.10| ), for all values of g and k. In this way, we find that the constants 
c 2 , c 3 and c 7 are uniquely determined as: 

c 2 = --; C3 = l; c 7 = -. (3.14) 

Note that the functions F 3 and F 7 , as well as the part of F 2 which determine the T 4 con- 
tributions [cf. Eq.( |2.5| )] are now uniquely fixed. This is in accordance with the argument 
that all hard thermal particles should contribute the same, up to a weight factor. The 
above relations imply further the equation: 
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14 



2k 2 F u + 



k 2 (k ■ q) 



(3.15) 



k 4 - 4 (k ■ q y 

Using ( |3.15| ), we see from Eqs. (|3.10| ) that the only independent functions left over are Fu 
and Fi 2 . From the general properties of the forward scattering amplitude, these functions 
must have the structure: 



Fn = en 



1 



+ 



k 2 + 2k-q k 2 - 2k ■ q 



(3.16) 



Fn — Cl2 



1 



(3.17) 



k 2 + 2k-q k 2 -2k-q)' 

where en and cu are constants which depend specifically on the nature of the thermal 
particles. For instance, in the scalar case we get: 



Cl1 = 12 ; 



° 12 ~ 36' 



(3.18) 



whereas in the case of internal gauge fields we find that: 



Cll = C12 = 0. 



(3.19) 



The above relations explain the features of the forward scattering amplitudes described by 
Eq. O [at f = -1/6 and m = 0] and by Eq. (|3~5|). 

IV. EXACT EVALUATION OF RADIATION FIELDS CONTRIBUTIONS 



We will now evaluate all finite-temperature contributions in closed form, using the tech- 



niques described in the first paper of reference [12] . To this end, we express the 1PI graviton 
2-point function in terms of the tensor basis ' Q/3 (u, K) in a way analogous to ( |2.9| ): 
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IF"- Q/3 (k) = J2 u i ^ K ) Tp al3 (u, K) . 



(4.1) 



i=l 



According to the discussion of the last section [cf. Eq. ( [2.101 )1, we can write the functions 
ilj (r, K) as follows: 
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IT (r, K) = k 2 



tt 2 T 4 
~T5 



-k (r,Ar) + |fc| 2 r^(r,Ar) 



(4.2) 



where the functions U (r, K) are given in appendix ([B]). Our task is to determine the func- 
tions Ni (r,K), which should be non-leading in the high temperature limit. For this, it is 
convenient to consider first the projections of the graviton 2-point function into the tensor 
basis 7T' Q/3 : 

Pi (r, K) = in""' aP (k) T. m a p{u,K). (4.3) 
Once we find these (see next), the functions IT in (|4.2|) can be determined by the relation: 



IT (r,K) 



liu, a/3 



T 



fiu, a/3 ) 



(4.4) 



where (T^-) -1 denotes the inverse of the matrix = T- lu ' al3 T j ^ a ^. 

We now proceed with the evaluation of the functions Pj (r, K) in [4.3| . From Table |, 
we see that for j = 4, 5, •••,14 these involve the contraction of Il^ lJ ' al3 (k) with i]^, r] a ^ 
or with the external momenta. Using the Ward ( [2.11| ) and Weyl ( j2.15Q identities, the 
corresponding functions Pj (r, K) will be given by a linear combination of graviton 1-point 
functions. These are proportional to T 4 [cf. Eq. ( [2.12| ) with m = 0], and so will contribute 
only to the functions li(r,K) in (|4.2| ). The functions Ni(r,K) are determined from the 
contributions corresponding to Pj(r,K) (j = 1,2,3), which are not proportional to T 4 . 
These contributions, which we denote by Pj (j = 1,2,3) can be found from Eq. ( |4.3|) by 
using for the graviton 2-point function the expression ( |1 . 1|) . Substituting here the expression 
for the forward scattering amplitude, we are lead to integrals of the form: 

1 r d 3 q Q s (I 1 



Is(k,T) 



+ 



(4.5) 



(2tt) 3 J 2Qexp(Q/T)-l\k 2 + 2k-q ' k 2 -2k- q / 

where S — 0, 2, 4 and Q = \q\. In terms of x = cos (9), where 9 is the angle between k and 
q, we find that the above expression becomes: 

1 k 2 r 1 dx r°° ,„ Q s+1 1 



(2tt) 2 4 









(k - 


k 


XJ 



dQ 



exp (Q/T)-lQ* + (27rTy) 



2 • 



(4.6) 
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where, 



1 



A' 2 



y 



AmT ko - 

Apart from simple functions, the integration in 



x 



(4.7) 



can be reduced to the basic integral 



i{y) 



QdQ 



1 



o Qi + (2TrT y yexp(Q/T) 



d 



£e[Re(y)] flny-i 



2y 



^{y))+ {y 



-y), 

(4.8) 



where ip(y) = — lnT (y) denotes the Euler psi function. The real-time limit of the Green's 
ay 

function can be obtained from the analytically continued imaginary-time formalism via the 
prescription k = (1 +ie)K , where e — > + and i^" is rea l- With the presence of the is 
factor being understood, we find in (|4.8| ) that Re(y) = e'Re(A;o), with e' — > + . 

Many of the angular integrations in ( |4.6|) can be easily done in terms of elementary 
functions, after changing variables from x to y. The most difficult one involve an integrand 
containing if) (y) multiplied by a power of y 11 , for n — 0, 1, 2, 3, 4. The relevant integrals can 
be put in the form: 



J n = e[Re(fc )] [Jn(*(fco)) - J„(-t(-A:o))] + [ko 



-k ] 



(4.9) 



where 



J n {t) 



' ■ m\ n-1 . 



My) ~tt- ^(y) 
2y 



dy. 



(4.10) 



The choice of C is immaterial, since any constant is irrelevant for our purposes because 



it cancels out in the expression (O). Here t (ko) and —t(—ko) denote the limits of the 



y-integration corresponding respectively to x = 1 and x = — 1 in Eq. (O). Hence: 



t(fc ) 



(4.11) 



These integrals can be expressed in closed form in terms of derivatives of generalized 
Riemann zeta functions ( {—n, t) for natural values of n = 0, 1, 2, 3, 4. For instance we have 
that: 
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Jo{t) 



vkT 



2y 



dy 



vkT 



tint - C' (0,t) - t- -hit 



(4.12) 



Ji(t) 



2/ 



dy = tC\0, t) - C'(-l, t) + ^ lnt - ^t 2 . (4.13) 



In the above expressions, the derivative is taken with respect to the first argument of the 
generalized zeta function. The functions J n (t) are discussed in more generality in Appendix 
C [cf. Eq. (|C6D 1. In this way, we find that the functions Pj (j = 1,2,3) in Eq. Q4.3|) are 
related to J n in (|4.9[) as follows: 



K 2 K 4 
Pi (r, K) = J , 



(4.14) 



„ K 2 T , s K 4 T K 2 K 2 T 

Pi = L (r) J n — r Ji Jo, 

2 48 V ; 32 4 1 2 2 ' 



(4.15) 



P 3 = 1 L (r) J — r Ji i-^ J 2 — r J 3 — J 4 , 

576 192 v ; 256 16 8 



(4.16) 



where we have defined: 



r r 4- 1 
L( r = - In— — - 1. 

y J 2 r- 1 



(4.17) 



With help of these relations, the functions iVj (r, if) can be explicitly determined from Eqs. 
Q4.2| ) and ( |4.4j ). After a straightforward calculation we obtain that: 



Aq(r, if) = Pi + if 2 P 2 + if 4 P 3 
N 2 (r, K) = K 2 P 1 + 2K 4 P 2 + 5K 6 P 3 
N 3 (r, K) = K 4 P 1 + 5K 6 P 2 + 35K S P 3 
N 4 (r, K) = -P x - K 2 P 2 + K 4 P 3 
N 5 (r, K) = -K 2 P 1 - K 4 P 2 + 5K 6 P 3 
N 6 (r, K) = -r (P x + 2K 2 P 2 + 5K 4 P 3 
N 7 (r, K) = -r (K 2 P 1 + 5if 4 P 2 + 35K 6 P 3 



(4.18a) 
(4.18b) 
(4.18c) 
(4.18d) 
(4.18e) 
(4.18f) 
(4.18g) 
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N 8 (r, K) = P 1 + (l + 2K 2 ) P 2 + (AK 2 + 5K 4 ) P 3 (4.18h) 
N 9 {r, K) = (2 + K 2 ) P 1 + (6K 2 + 5K 4 ) P 2 + (30K 4 + 35K 6 ) P 3 (4.18i) 



iVi (r, K) = K 2 P 1 + [3K 2 + 5K 4 ) P 2 + [30K 4 + 35K b ) P 3 (4.18j) 

N u (r, K) = -r [P x + (2 + 5K 2 ) P 2 + (20K 2 + 35K 4 ) P 3 ] (4.18k) 

N 12 (r, K) = P 1 + (A + 5K 2 ) P 2 + (8 + 40K 2 + 35K 4 ) P 3 (4.181) 

N 13 (r, K) = -Pi - K 2 P 2 + (4iT 2 + 5K 4 ) P 3 (4.18m) 



JV 14 (r, X) = r (P + iT 2 P 2 + 5K 4 P 3 ) (4.18n) 

At high temperatures, the Riemann zeta functions can be expanded in a power series 
in t, as shown in Eq. (|C10| ). Then, in the high temperature domain, we can express the 
functions Pj (j = 1, 2, 3) as a series of powers of {1/T). The dominant terms in these series 
are given by: 

K 2 

P^-j 6 (4.19) 
Pi = |g £M (4.20) 

3 576 y 1 

Although the above contributions are gauge invariant, they do not directly describe the 
physical properties of the plasma at finite temperatures in a gravitational field. This problem 
is related to the fact that the thermal graviton 2-point function depends on the choice of 
the basic graviton fields. 

V. THE GRAVITON SELF-ENERGY AT FINITE TEMPERATURE 

In thermal quantum field theory the 1PI contributions to the graviton 2-point function 
are in general dependent on the parametrization of the graviton fields. However, as shown 
in the second work of reference ||12||, the traceless quantity: 



u ^,a(3 ^ = n ^,a/3 _ ]_ (^^P ^ + p ™ + rf a W p ^ + T]^U P p ^ a ) , (5.1) 
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represents at high temperatures a quantity which does not depend on the choice of basic 
graviton fields. In this domain, the masses can be effectively neglected so the theory is 
invariant under conformal transformations. The representation-independence of IT then fol- 
lows in consequence of the Ward and Weyl identities. As we have seen, the contributions 
from internal massless particles, are invariant under local coordinate and conformal trans- 
formations. Consequently, the physical amplitude given by Q5.1|) can be identified in this 
case with the graviton self-energy even at finite temperature. 

However for contributions from thermal matter, which are characterized by the presence 
of massive particles, IP 1 '' Q/3 is no longer independent of the graviton field parametrization 
at finite temperatures. Our task is to generalize (|5.1|) in such a way that the corresponding 
quantity should represent a physical graviton self-energy at all temperatures. To achieve 
this we consider the effective action which generates the one-particle irreducible thermal 
Green functions. In the representation ( |2.2j ) for the h fields, we have that: 



S eff = Vfh ap (0) + ~ f d 4 Hir a/? (k) V (*) K P (-k) + ■■■. (5.2) 

Starting from these fields, the most general re-parametrization of the graviton fields can be 
written as: 

= ahT + b h x x jf v + c h\ h" u + d h\ + e (h\) %^ + / h a ph a(3 if + • • • , (5.3) 



where a, b, c, d, e, and / denote arbitrary constants. For example, a basic graviton field 
often used in the literature which is defined by [|] 



y/-9 W v = ^ v + i&r, (5.4) 
corresponds to a special case of (|5.3|), with: 

a = — 1; b = -; c = — ; d = 1\ e = -; f — ——. (5.5) 
> 2 2 8 4 y 1 

In what follows we shall assume, without loss of generality that a 2 = 1. This can always be 

achieved by a further rescaling of the h, fields in ( |5.3j ) [see ref. Furthermore, we shall 



consider for simplicity the class of parametrizations characterized by the conditions 
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ab-ac + 2bd = 0; b + 2 f = 0, (5.6) 
which are explicitly verified by all the graviton representations discussed in the literature 



Since the effective action is invariant under a general re-parametrization of the graviton 
fields, it can be written in terms of h, as: 

S eff = rtfh laf) (0) + \f d A kU^ a/3 (k) V (k) h, a/3 {-k) + ■■■. (5.7) 



Identifying with the help of (|5.3| ), the corresponding terms in (|5.2|) and (5/7 ), we obtain the 
following relations: 

ir = oir+6iri%p, (5.8) 

nr aP (k) = K; a "(k) +ab (nc p , rf* + ng> %^) + b 2 Kfx u ^ 

*/3 iuA 



+ c(rc^ + r^ ; V 
d 



(5.9) 



+ 2e T p h/ fl^rf? + fV p hip {rj^r] va + rfrf p \ 

When the theory is invariant under conformal transformations, the following Weyl iden- 
tity holds pf: 

jr P ,nv _ a + ^c + 2d _ (k^q) 
U hlP ~ a (cl -\- 46) Lh ' ~ (5 ' 1Uj 

where is a traceless function. In general, this is no longer true in the presence of thermal 
matter at finite temperature. In order to take this fact into account, we generalize (|5.1f) by 
considering the following traceless quantity: 



(5.11) 



^rj fJ,a A u/3 + r]^ 13 A va + rf a A^ + rj^A^ 01 ^ 
where the tensor A Miy is given by: 

br = \ (if/* + f") - i (if/ + f" p ) ir ■ (5.i2) 
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We remark that when the Weyl identity (|5.10|) is applicable, A MW vanishes so that Eq. (|5.11|) 



reduces to (|5.1|) as expected. For this reason, only the contributions from thermal matter 
will appear in Eq. fl5.12| ). 

It is now straightforward to verify, with the help of the relations ( |5.6| ), ( |5.8| ) and ( p.9| ) 
that: 

n{f'^(jfe) = n%' a/} (k). (5.13) 



This equation shows that the graviton self-energy given by the relations (|5.11|) and ( 5.12Q 



is invariant under re-parametrizations of graviton fields at all temperatures. In order to 
understand the mechanism which enforces the above property, using the relations ( |5.1D , 
( p. 11| ) and (|5.12 ), we write the expression for the graviton self-energy in the form: 



(jfc) = n^- Q/3 (fc) + n^ Q/3 (k) , (5.14) 

where 11^ is given by: 

n^ a/3 (k) = - {rfv v ^ + n^f va + r] ua r^ + 7f£i>*) 

As mentioned before, A vanishes in the case when the thermal fields can be considered as 
being effectively massless. The contributions associated with the graviton 1-point function 
in ( |5.15| ) can be represented diagrammatically as shown in Fig. |5|. Both II and T depend 
individually on the choice of basic graviton fields in a way that ensures IT to be independent 
of these parametrizations. Hence, in order to obtain a physical self-energy, one must consider 
in addition to 1PI graviton 2-point function, also the corresponding "tadpole" contributions. 
Since the graviton self-energy (|5.11|) is parametrization-independent, it may be conveniently 
evaluated in the representation (|2.2|) of the gravitational fields, from the contributions of 
thermal matter and radiation fields given respectively by Eqs. (|2.9|) and (|4.1|). 
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VI. THE EFFECTIVE GRAVITON PROPAGATOR 



In order to investigate the thermal mass of gravitons, we will study the properties of 
the poles in the effective graviton propagator. This is obtained by iterative insertions of the 
physical self-energy in the classical graviton propagator k~ 2 P^, where: 



r aP — 2 = a P \°- L ) 



which is insensitive to changes of parametrizations ||[12 



As we have seen, because of the inclusion of graviton 1-point functions, the graviton self- 
energy ( |5.14j ) is also independent of the parametrizations of graviton fields. These properties 
ensure that physical quantities such as masses are independent of the choice of basic graviton 
fields. Using the fact that the physical self-energy is traceless, P behaves effectively like the 
identity when acting on IT. Hence, the effective graviton propagator can be written in the 
form: 



i ~ / 1 \ 2 

' ^ ' (6-2) 



PS + ^ {k) + ) WW) + • • • 



The right-hand side of this equation sums up to a geometric series, giving the relation: 

(tfpg _ ft^) d% = (6.3) 

The effective propagator satisfies certain fundamental constrains. In view of the traceless 
property of IT, Eq. (|6.2|) requires that: 

D p p '" v = (6.4) 

Furthermore, the Ward identity ( p.llj ) expressing gauge invariance requires a longitudinal 
contribution in II connected with the background energy- momentum tensor [cf. Eq. (2.16) 
in the second paper of Ref. ||12|| 1. Considering for definiteness the high-temperature limit, it 
is then straightforward to verify that Eq. (|6.2|) implies: 

= ^ - \^ + - 4 «)(^ 2 1 - ti)-l, v (6.5) 
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where the energy density p is given by: 



P 



CU7T 2 T 4 

30 



(6.6) 



Here to denotes the total number of degrees of freedom of the thermal particles. 

In what follows we shall be interested only in determining the effective dynamical masses 
in the static case ko = 0, which are relevant in the process of dynamical screening. To this 
end, we project the corresponding contributions of II into the following traceless normal 
modes: 



rpfllS, afi 



rp[if, a/3 



1 



1. 



-3 W - -Uu) + -T 5 (u) 



/xu, a(3 



--T 2 (u) + 2T 3 (u) 



liv, a(3 



(6.7) 
(6.8) 



where the tensors T^"' (u) (i = l,---,5) are obtained from the corresponding ones in 
Table | by replacing X a with u a = 5$. The normal modes Tj and T$ are idempotent (up 
to a minus sign) and orthogonal to each other. While the mode Tj is three-dimensionally 
longitudinal, the mode T$ is spatially-transverse in the sense that: 



h h-T 13 ' ^ v — n- 

ru % hjjJ. g — U, 



(z, j = 1,2,3) . 



(6.9) 



In terms of these tensors, we can decompose the self-energy as follows: 



n""^(jfco = 0) = k 2 p (rp^ - brp<*p 



(6.10) 



It is then easy to invert Eq. ( |6.3|) , yielding for the effective graviton propagator the result: 



1 



1 



-T, 



f)^,*? r k = 0) = T y,ap _ ; 

k 2 k 2 I p ' 

where the normal mode Tjj given by: 



fiv, a/3 



+ 



k 2 



-T 



J 



T, 



U 



-~T x {u) + -T 2 (u) - -T 3 (u) + -T A (u) - -T 5 (u) 



/iv, afi 



(6.11) 



(6.12) 



is orthogonal to the modes Tj and T5. 
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We see that in the normal mode Tu, the gravitational plasma is unscreened. This is 
somewhat similar to the spatially-transverse mode in the QCD plasma. On the other hand, 
a non- vanishing screening mass appears in the mode T5: 

ml = ^P = ^l (6.13) 
3 3 

analogously to the behavior shown by the spatially-longitudinal mode in the QCD plasma. 
The mode Tj is characterized by an imaginary mass: 

m] = -32nGp (6.14) 

which is similar to the classical Jeans mass. This anti-screening mode indicates a gravita- 
tional instability for density fluctuations with wavelength larger than |mj| _1 , owing to the 
attractive nature of gravity. One may generalize this calculation by including internal gravi- 
tons in thermal equilibrium at high temperatures 0. Their contributions will not affect the 
above conclusion, since these change only the weight factor w appearing in p [cf. Eq. (|6.6| )1 
which counts the total number of degrees of freedom. 

In conclusion we emphasize that this work, like that of Ref . || , has been concerned with 
gravitational perturbations at finite temperatures around the Minkowski background, in an 
asymptotically flat space. Using a different approach, based on the study of small distur- 
bances around the solutions of Einstein equations, Rebhan || has performed a rather com- 
plete investigation of the gravitational instabilities. These metric perturbations are relevant 
at high temperatures in the context of a radiation dominated Robertson- Walker universe. 
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APPENDIX A: 



In this appendix we present the Feynman rules for the couplings and propagators in- 
volving scalar, gauge and graviton fields. These rules can be obtained from the respective 
Lagrangians given in Eq. fl2.1|) and Eqs. Q3.1J) , (|3.2| ) and (|3~3"1). In a perturbative calculation 



we first have to expand all the metric dependent quantities up to some given order in the 
graviton field h. These expansions and the subsequent reading of the momentum space 
Feynman rules is a straightforward procedure (but a very tedious task for humans) which 
was accomplished using a algebraic computer algorithm written in Mathematica. Here we 
will only present the results for the vertices involving up to two gravitons, which are rele- 
vant for the calculation of the graviton polarization tensor. We will also restrict only to the 
Abelian couplings of the gauge fields [cf. Eq. Q3.1J) ]. 

In all the expressions which follows we will always denote the graviton momenta and 
indices by [kx, (/x, u)] and [k 2 , a, /?)]• The momenta of scalars and ghosts are denoted by 
Pi and pi- The gluon momenta and indices are denoted by \p%, p] and [p 2 , o\. Using this 
notation, we obtain from the Lagrangian (|2.1| ) the scalar- scalar- graviton interaction vertex 

-V^ calar (ki;p 1 ,p 2 ) = P w P 2 u+ P 1 uP 2f ,-pi-P2V^-m 2 r] IIiy + 2^ (k lfl k lV - fc?^) (Al) 

and the scalar- scalar- graviton- graviton interaction vertex 
16 

— V^ffiih, k 2 ;pi,p 2 ) = -8p lU p 2 pr] a ^ + 2p 1 -p 2 r] afl r] (3u + 4:p ia p 2 pr]^- 

Pl ■ P2 Va(3 Vfiu + m 2 (2 Vafi Vpv ~ Vafi 7]^) + 

f (4 k lfl k lU 7] a p + 4 k x fj, k 2 y 

6 k 1 p k 2V r] afl — 4 k 1 @ k lfl r] au + 2 k 1 p k 2fl n au — 

(A2) 

4 k ia k lU 7]^ - 2 k ia k 2U n^ + 2k x -k 2 n ai/ r)^- 
4 k ia k 1 ^ l rjp u 8 k 1 ^ l k 2Ct f]p v + 2 k lCt k 2fM T)p u -\- 
8 k\ r] atl T] Pv +Ak x -k 2 r] ail rip u + 8 k ia k l/3 r]^+ 
4 krf k 2a T]^ -Akf r] aP -q^ -2k x -k 2 r] al3 r]^ u ) , 
where in the expression above one has to perform a symmetrization over the graviton indices 

and permutation of the scalar particles. 
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From the ghost Lagrangian we obtain the ghost- ghost- graviton interaction vertex 
2 

-V^ h ° st (h; PUP2) = ~PiuP 2 p-p 1 pP 2 u+Pi -P2V^, (A3) 

and the ghost- ghost- graviton- graviton interaction vertex 

8 2 

~ v ,J%p( k i, h] PhP%) = -4p lM p 2i , Va/3 + 8p lP p 2U 7] Q ^ -2p 1 -p 2 r] afl rj Pu +p 1 ■ p 2 r] a/3 r]^. 

(A4) 

Expression ( |A4j ) has to be symmetrized over the graviton indices. Notice that, as can be 
easily seen from Eqs. (|2.1j) and ( p.3|) , the interaction vertices of ghosts or scalars with the 
graviton field, differ only by a minus sign when £ = m = 0. The corresponding propagators 
are given by 

!1 
— for the ghost, 

k\ . (A5) 

— — for the scalar 
m z — k 2 

From Eqs. (|3.1| ) and (|3.2j ) we can obtain the gauge fields Feynman rules in the general 
covariant gauge characterized by the gauge fixing parameter a. In this class of gauges the 
gauge field propagator is given by 



k k 

Vnu -(1-0) -^r 



(A6) 



The interaction vertices will also depend on the parameter a. The gauge- gauge- graviton 
coupling is 

4 1 

- V ^ 9 a( k i;Pl,P2) = -PjaPipV^u + P 1 aP 2 uV^ + PiuP 2 pV^ + 
K 

Pxa P 2 p Vup + P w P 2 p Vua ~ 2 Pi ■ P2 V«tr~ 

2 P 1 p P 2 v Vpa +P1-P2 Vpu Vpa+ ( A7 ) 
- (-ZPipPiaVpv ~P lP P 2 aVpu ~ 2PtvP 2 *Vw+ 
^P!uP iP Vpa + 2PipP 2 vVtlv) 7 

and the gauge- gauge- graviton- graviton vertex is 
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— V ^a/3-p,a( k l ,h]Pl,P2)= 2p ic p 2f) T] ap T]p„ - 4 p ia p 2V T] ap T)p p - 4 p ±v p 2p T] afl T]p a + 

/ tPinP2vVapV(3cr -p i aP 2P Vaf3Vpu + 2p ltT P a pVapVl^ + 
2 Pi/9 P 2 p Vaa Vpu + 2 Pi <J P 2 a V/Sp V/JV + 2 Pi* P 2 p Vficr ?/> _ 

4 Pi ■ p 2 ^ - 4 p 1Q jv 77^ 77 w - 4 p 1<r p 2/3 77^ 77^- 

4 Pi/3 P 2A i ^a<x ?7i/p + 8 pi ■ p 2 T) afJl Tjfa T) up - 4 p i/3 p 2p 7] ap T] ua + 

4p ia p 2 /3 Vpp Vua + 8p lU p 2 p r] ap r] pa -2p 1 -p 2 r] ap rjp v r] pa - 

4 P ia P 2 (3 Vpu Vpa +P1-P2 Va(3 Vpu Vpa+ 

- (16 k lU p ip r] aa rjpn - 8 k ia p ip r] ap <qp v - 2 p ip p 2a r] ap 77^+ 
2 k lP k 2a rj a p rj pu + 4:k 1(T p lP rj a pri pu + p lP p 2(T r? Q/3 77^- 
8 P^ ^ap ?/> - 4p 1( g p 2a 7] ap i] pv - 4 k lP k 2/3 rj au i] pv - 

8 p ip T] aa Tjfja, - 4 p ip p 2/3 77 a(T 7/^ -Ak lU k 2U 77 a/3 7/ MP + 
8 Pi/3 P 2 a Vau Tj pp + 8 k lV k 2 (3 T] aa T] pp + 8 p iu p 2 p T] aa T] pp - 

8k lU p lP T] a/3 rjfur + 8 fc 1/3 p lP 77^ 77^ + 8 p lP p 2/3 r7 Q! , 77^+ 
16 k xV p lf 3 ri a p rj pa + 8 k ia p lP r)p v rj pu ) . 

(A8) 

Similarly as in the scalar vertices, one has to symmetrize the gauge field vertices over the 
graviton indices and include the permutations of the gluons. In all interaction vertices there 
is momentum conservation, with all momenta inwards. 

From the expressions above we can perform the explicit computation of the scattering 
amplitudes shown in the figures |2] and f|. This was done using these Feynman rules as an 
input to a algebraic computer program. 

APPENDIX B: 

In this appendix we present the leading and next-to-leading structure functions for the 
matter contribution to the polarization tensor. The leading structure functions presented 
here are the same for all thermal particles. The explicit result for the functions li(r,K), 
rii(r,K) and Si(r,K) appearing in Eq. (|2.10|) is the following: 

24 



= -5K 4 5K 2 £ , 5K 6 L 

1 192 16 "r" 64 

c — _5 _L 5X 2 L 

S 1 ~~ 8 ^ 16 



2 3 + 12 24 8 

_ -25 K 6 _,_ 5if 6 L _,_ 25K S L 
2 192 ^ 32 ' 64 

_ 5 _ _,_ 15if 4 L 

2 8 16 ^ 16 



r -K 2 _|_ 7 if 4 _ 35 K 6 _|_ 35 K 8 L 

3 3 12 24 ^ 8 

-5K 6 _ 175 K 8 _| 25 X s L , 175 if 1 " L 
3 32 192 ^ 16 64 

_ 5 if 2 _ 25 K 4 _,_ 75 K 6 L 
3 8 16 t 16 



—K 2 , K 4 L 
24 8 

—5 if 2 _ 5^ _ 5X 2 g , 5 K 4 L , 5if 6 L , 5K 4 £L 

32 192 8 "l" 16 T 64 T 4 

5if 2 L 
16 



if 2 _ 5 if 4 , 5if 6 L 
12 24 8 

-5K 4 _ 25 if 6 _ 5if 4 g , 5if 6 L i 25 K 8 L , 15-ff 6 gL 
32 192 8 ~ l ~ 8 64 8 

-5 if 2 _,_ 15if 4 L 
16 16 



7 - f =1 _L SiC 2 . _ 5K 4 L \ 
6 \ 12 24 8 J 1 

( 25K 4 5K 4 L 25K 6 L \ 
' i6 V 192 32 64 J 

66 V 16 16 J ' 



j fl _ 7K 2 _,_ 35 K 4 _ 35 K G L 

7 V 3 12 24 8 

— ( 5I<4 _|_ 175 if 6 _ 25 K fi L _ 175 K 8 L \ 
7 V 32 192 16 64 J 

o _ ("-5 I 25if 2 75if 4 L \ „ 
6 7 - V 8 16 16 j 7 
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; __L 5K 2 i K 2 L i 5K 4 L 

8 12 24 2 8 

-5K 2 _ 25 -ftT 4 _|_ 5_£ _,_ 5.R: 2 L i 15-ft~ 4 L i 25X 6 L 
8 48 192 + 16 + 32 32 64 

— _JL _|_ 5L _|_ 15X 2 L 
8 16 ^ 8 16 



_ 1 _ 2X 2 _ 35 if 4 , 15X 4 L | 35 K 6 L 
9 6 3 24 "I" 4 8 

— 15 -ft" 4 _ 175 -ft" 6 , 5K 2 £ _,_ 15-ftT 4 L , 125 if 6 L , 175 K 8 L _ 15 K 4 g L 
-"9 16 192 8 16 32 64 8 

— _5 _ 25 -ft -2 , 15-ftT 2 -L i 75 K 4 L 
9 8 16 4 16 



/ — I _ 2 -ft -2 _ 35 -ft" 4 _|_ 15-ftT 4 L _|_ 35 K 6 L 
1 10 — 6 3 24 ~r 4 ~T 8 

-5X 2 _ 15 -ftT 4 _ 175 -ftT 6 , 45 K 4 L , 125 K 6 L , 175 K 8 L 
' L 10 32 16 192 32 32 64 

_ _5 _ 25 -ft -2 , 15X 2 L _,_ 75 K 4 L 
b 10 ~~ 8 16 ~r 4 16 



(l , 35 X 2 _ 5 -ft" 2 L _ 35 X 4 L \ 
U 24 2 8 / 

( 65 -ft -2 , 175 X 4 _ H _ 5K 2 L 
\ 96 192 8 8 

/ 25 _ 15L _ 75 K 2 L \ 
V 16 8 16 J 



25 K 4 L _ 175 -ft" 6 L , 15-ft: 2 gL 



Z 12 = -g - M! + L + 5 ^ i+ 3i|lL 



_ __5 115 K 2 _ 175 -ft" 4 _i_ M _|_ 15X 2 L , 75 X 4 L , 175 X 6 L _ 5gL _ 15.ft: 2 gL 

12 24 96 192 f 4 + 8 + 16 f 64 2 4 

_ _25 _ _5 I 15L _|_ 75 K 2 L 

b 12 ~~ 16 8if 2 + 4 + 16 



_ _ J 5K 2 , X 2 L , 5K 4 L 

13 12 24 2 8 

-25 X 2 _ 25 X 4 _ 5X 2 g , 5K 2 L , 15X 4 L , 25X 6 L , 5 K 2 g L . 15 K 4 £ L 
' 13 96 192 8 + 16 + 16 + 64 + 8 + 8 

— _JL _|_ 5_L , 15.ft: 2 L 
13 16 + 8 + 16 



7 _ /-l i 5X 2 5X 4 lA 
*14 ~ U2 + ~24 S - J r 

/ 5_ft^ i 25 K 4 , 5X 2 g _ 5-ft~ 4 L _ 25 K 6 L _ 15X 4 gL 
II 14 ^ 32 192 8 8 64 8 

s 14 _ \ie 16 ) r 
The dimensionless quantity L is a function of r given by: 
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APPENDIX C: 



Here we calculate the integrals 

Jn(t) = 



. rr \ 71—1 
11X1 \ rt 



\k\ 



dyy n 



2y 



(CI) 



in terms of the generalized ( function, defined as 

oo ^ 

c(2iS)= S(i^ 



(C2) 



To this end we express the function as 

1 



£-+0 



C(l + e,2/) 



(C3) 



and use the formula 



dy((z,y) 



1 



l-z 



C(l-z,t), 



(C4) 



that can be easily verified from eq. (|C2| ), and can be generalized to 

A(n-Z)! r(l-z) 



J dyy n ((z,y) = J2 



-t n ((z-l-l,t), 



(C5) 



^ n\ Y(l + 2-z) 

which is eq. ( |U4|) integrated by parts n times. 

Substituting eq. ( |C5D in eq. ( [CI] ) and using eq. QC5| ) we can verify, using the properties 
of ( function, that the divergent term as e — > cancels out, as expected. For n 7^ the 
remaining terms give: 

t n+l 



Jnif) 



1*1 / 



n + 1 



^-i-s<-^(;ifw') + 



E^f")fti|^ +1 ( f )}, 



(C6) 



where B n are the Bernoulli polynomials ||11|| . For n — 0, Jo(£) is given by Eq. ( |4.12| ). 

Now we discuss the behavior of the generalized zeta function for asymptotic values of the 
parameter t(ko) = i k °^^ , which correspond to high temperature expansion. To this end, 
we start from the representation [11] : 
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COM) 



l l 



t z T(z) Jo e 



CO r£ Z lg 



fix. 



(C7) 



Expanding flC7|) in power series of t, making use of the integral representation of Rie- 
mann's zeta function and Euler's gamma function we find 



a/ 1 ™r(z + l)(-t)' 



* 2 S r(*) zi 

Taking the derivative of flCSj) with respect to 2; we obtain in terms of the psi function 
ip(z) that 



CCM) = £ 



H)T(z + o 



{[^ + /) _ ^)] C(z + /) + C \z + Z)} - t*hx{z). (C9) 



We are actually interested in the values of ('(z,t) for z — > — n where n is a natural 
number. After a straightforward calculation we obtain 



•n 1 

fn+l 



C(/-n)-C(i-n) £ 



k=n-l+l 



k 



t l - t n \n(z) 



n + 



\ /l — 1 / / — TL~ \~ 2 



(CIO) 



With help of this formula we can compute the functions J n in Eq. ( |C6|) and express 
P n from Eqs. (|4.14 ), (|4.15|) and ( 4.161) as a series of decreasing powers of T. Then, it is 
straightforward to arrive at Eqs. ( 4.19 ), ( 4.20 ) and (|4.21 ). 
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TABLES 



Tf aP (X, Y) = rf v rfv + n a * *f v 

T iw, = xv (X? rf v + X a 77^) + X v (x$ rfv + X a r/^) 

T p Y ^ = X a X {1 xfl X u 

T£ u ' aP (X,Y) = rf^rf v 

T txv, Y } = X v X v T] a @ + X a X? rf v 

Q/3 (A, Y) = X? (Y u rfv + YV rf v ) + Y p (X u rfv + XV rf v ) 
+ X a {y v rfv + YV rf v ^ + Y a (x v rfv + XV rf v ^ 
Tj U ' a(3 (X, Y) = Y u X a X? XV + YV X a X$ X v + X a XV X u + Y a X? XV X 
rpgv, a P^x, Y) = Y 13 Y v rfv + Y 13 YV r] au + Y a Y v rfv + Y a YV rf v 
Tp aP {X, Y) = YV Y v X a X? + Y a Y$ XV X v 
T$' Q/3 (A, Y) = (y? X a + Y a X^ (Y u XV + YV X v ) 

rpixv, ^Pt x ^ Y) = Y@ YV Y u X a + Y a YV Y v X^ + Y a Y@ Y u XV + Y a Y@ YV X u 

T} xv, a/3r x ^ Y ) = Y a Y 13 YV Y v 

Ti xv, a p^ Y } = Y v y v rfP + Y a Y? if 

T} xv, af3r x ^ Y } = (yv Xii + Y v x u ) rfP + {y? X a + Y a X^ rf 



TABLE I. A basis of 14 independent tensors Tf' ap (X, Y) 
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FIGURES 

FIG. 1. Lowest order matter contributions to the thermal 1PI graviton two-point function. 
Curly lines denote the external gravitational field and solid lines represent the scalar particle. 

FIG. 2. The forward scattering graphs corresponding to Fig. 1. Crossed graphs with (k <-> —k) 
are to be understood. 

FIG. 3. One-loop contributions of radiation fields to the graviton polarization tensor. Wavy 
lines denote the gauge field and broken lines represent ghost particles. 

FIG. 4. Forward scattering diagrams containing ghost particles connected with Fig. 3. Crossed 
diagrams (&•<-> —k) should be included. 

FIG. 5. Lowest order contributions of graviton 1-point functions to litad- The black dot repre- 
sents terms proportional to rj. 
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